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Floppy modes and non-afRne deformations in random fiber networks 
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We study the elasticity of random fiber networks. Starting from a microscopic picture of the non- 
afpne deformation fields we calculate the macroscopic elastic moduli both in a scaling theory and 
a self-consistent effective medium theory. By relating non-afSnity to the low-energy excitations of 
the network ( "floppy- modes" ) we achieve a detailed characterization of the non-affine deformations 
present in fibrous networks. 
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Materials as different as granular matter, colloidal sus- 
pensions or lithospheric block systems share the com- 
mon property that they may exist in a highly fragile 
state [l|, y| ■ While in principle able to withstand static 
shear stresses, small changes in the loading conditions 
may lead to large scale structural rearrangements or even 
to the complete fluidization of the material [1,1^^. To 
understand the extraordinary mechanical properties of 
these systems new concepts have to be developed that go 
beyond the application of classical elasticity theory and 
that sufficiently reflect the presence of the microstruc- 
ture Q. One example is the "stress-only" approach to 
the elasticity of granular materials ^ , where the elimina- 
tion of the kinematic degrees of freedom accounts for the 
infinite stiffness of the grains. This seems to capture the 
inhomogeneous distribution of stresses in the sample and 
their concentration along the so called force-chains 0- 
In jammed systems of soft spheres, on the other hand, 
fragility has recently been shown to directly affect the 
deformation response of the system. While it may in- 
duce anomalous deformation fields that strongly deviate 
from the expectations of homogeneous elasticity ("non- 
affine" deformations) it may also lead to a prolifer- 
ation of low-frequency vibrational states far beyond the 
usual Debye-behaviour of ordinary solids |^ . It has been 
argued that these low-energy vibrations derive from a set 
of zero-frequency modes (floppy modes) that are present 
just below the jamming threshold and relate to the 
ability of the structure to internally rearrange without 
any change in its potential energy. This concept of floppy 
modes has also been used in connection with elastic per- 
colation networks where the fragile state is reached by 
diluting a certain fraction of nearest-neighbour contacts. 
In these systems constraint-counting arguments may be 
used to determine the percolation transition at which the 
system looses its rigidity [Tl) . 

Here our focus is on a particular class of heterogeneous 
networks composed of crosslinked fibers. These systems 
have recently been suggested as model systems for study- 
ing the mechanical properties of paper sheets Il2l or bio- 
logical networks of semifiexible polymers [l^l Il4l|. While 



these networks are known to have a rigidity percolation 
transition at low densities [isl H^ . we show here that 
even networks in the high-density regime in many ways 
resemble the behaviour of fragile matter, despite the fact 
that they are far away from the percolation threshold. 
We identify the relevant floppy modes and highlight their 
importance for understanding the macroscopic elasticity 
of the network. In particular, we will be able to ex- 
plain the occurence of an anomalous intermediate scaling 
regime observed in recent simulations 0, 0, . In this 
regime the shear modulus was found to depend on den- 
sity (measured relative to the percolation threshold) as 
G ^ 5p^ with a fractional exponent as large as ii « 6.67 
|0| . Also, highly non-affine deformations [tR 0| as well 
as inhomogeneous distribution of stresses in the network 
have been found. Heuristic non-affinity measures have 
been devised 0, 0| , however, little is known about the 
actual nature of the deformations present. While the 
expression "non-affine" is exclusively used to signal the 
absence of conventional homogeneous elasticity, scarce 
positive characterization of non-affine deformations has 
been achieved up to now . This Letter tries to fill this 
gap by characterizing in detail the non-affine deformation 
field present in fibrous networks. By relating non- affinity 
to the floppy modes of the structure we can, starting from 
a microscopic picture, calculate the macroscopic elastic 
moduli both in a scaling theory and a self-consistent effec- 
tive medium approximation. In analogy with the affine 
theory of rubber elasticity for flexible polymer gels, our 
approach might very well serve as a second paradigm to 
understand the elasticity of microstructured materials. 
Due to the proximity to the fragile state, it might also be 
of relevance to force transmission in granular media and 
to the phenomenon of jamming. 

The two-dimensional fiber network we consider is de- 
fined by randomly placing N elastic fibers of length If 
on a plane of area A — such that both position and 
orientation are uniformly distributed. We consider the 
fiber-fiber intersections to be perfectly rigid, but freely 
rotatable cross-links that do not allow for relative sliding 
of the filaments. The randomness entails a distribution 
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of angles 6*6 [0, tt] between two intersecting filaments 

= (1) 

while distances between neighbouring intersections, the 
segment lengths Ig , follow an exponential distribution [2ll | 

P{h)^{Q-'e-'^/^'^K (2) 

The mean segment length (Ig) is inversely related to the 
line density p = Nlf/A as (Ig) — 7r/2p. The segments 
are modeled as classical beams with cross-section ra- 
dius r and bending rigidity k. Loaded along their axis 
( "stretching" ) such slender rods have a rather high stiff- 
ness fc|| (Is) — 4K//sr^, while they are much softer with re- 
spect to transverse deformations k±{ls) = in/tl ("bend- 
ing"). Numerical simulations for the effective shear mod- 
ulus G of this network have identified a cross-over scaling 
scenario characterized by a length scale ^ = lf{6plf)~^ 
and v « 2.84 16] that mediates the transition between 
two drastically different elastic regimes. For fiber radius 
r 3> ^ the system is in an affine regime where the elastic 
response is mainly dominated by stretching deformations 
homogeneously distributed throughout the sample. The 
modulus in this regime is simply proportional to the typ- 
ical stretching stiffness, Gaff oc fc||((Zs)) and independent 
of the fiber length If. This is in marked contrast to the 
second regime at r <C ^. There, only non-affine bending 
deformations are excited and the modulus shows a strong 
dependence on fiber length Gna oc k±{{ls)){lf / {ls))'^~^ ■ 
Using renormalization-group language the parameters r 
and / f may be viewed as scaling fields (measured in units 
of the "lattice-constant" (Is))- The stretching dominated 
regime may then be characterized by an (afhne) fixed- 
point at — > oo and finite radius r ^ 0. On the other 
hand, the (non-affine) fixed-point of the bending domi- 
nated regime is obtained by first letting r —* and then 
performing If oo. This suggests that the elastic prop- 
erties in the latter regime may be analysed at vanishing 
radius r = 0, that is by putting the system on the stable 
manifold of the fixed point. 

In the following we will exploit this limit to calculate 
the modulus Gna in the non-affine regime. Central to the 
analysis is the recognition that in this limit the ratio of 
bending to stretching stiffness fcj_/fc|| oc tends to zero 
and bending deformations become increasingly soft. We 
thus obtain the much simpler problem of a central-force 
network. However, as only two fibers may intersect at a 
cross-link the coordination is z < 4 |2^ and rigid regions 
may not percolate through the system |23 . lisj . This 
implies that on a macroscopic level, the elastic moduli 
will be zero, while microscopically displacements can be 
chosen such that segment lengths need not be changed. 
These are the floppy modes of the structure that entail 
the fragility of the network in the bending dominated 




FIG. 1: Construction of a floppy mode by axial displacement 
Sz of the primary fiber (drawn horizontally) and subsequent 
transverse deflection y — ~Sz cot 6 of tlie crosslinks to restore 
the segment lengths on the secondary fibers (dashed lines, 
possible to first order in 5z). Initial cross-link positions are 
marked as black squares, final configurations as red circles. 



regime. It has been argued that a critical coordination 
of Zc = 4 is necessary to give the network rigidity |23|. 
This value defines the "isostatic" point, which in our net- 
work corresponds to taking the limit // — > oo. Thus, 
we arrive at the conclusion that isostaticity and the on- 
set of rigidity seem to be intimately connected to the 
fixed point governing the non-affine regime. While it is 
usually not possible to deduce the specific form of the 
fioppy modes, the fibrous architecture allows for their 
straightforward construction. In a first step (see FigQ 
we perform an arbitrary axial displacement Sz of a, given 
(primary) fiber as a whole. This, of course, will also af- 
fect the crossing (secondary) fibers such that the lengths 
of interconnecting segments change. In a second step, 
therefore, one has to account for the length constraints 
on these segments by introducing cross-link deflections 
yi = —Sz cot 9i transverse to the contour of the primary 
fiber. As a result all segment lengths remain unchanged 
to first order in Sz 26]. The construction is therefore 
suitable to describe the linear response properties of the 
network, while at the same time it offers an explanation 
for the stiffening behaviour found in fully nonlinear sim- 
ulations [l9l |24| . Any finite strain necessarily leads to 
the energetically more expensive stretching of bonds and 
therefore to an increase of the modulus. 

The identified modes take the form of localized exci- 
tations that affect only single filaments and their imme- 
diate surroundings. By superposition we may therefore 
construct a displacement field that allows the calcula- 
tion of macroscopic quantities like the elastic moduli. To 
achieve this we need to know the typical magnitude of 
displacements Sz of a given fiber relative to its surround- 
ings, the crossing secondary fibers. Since Sz is defined 
on the scale of the complete fiber we do not expect any 
dependence on average segment length (Ig), such that 
Sz (X If remains as the only conceivable possibility. Al- 
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ternatively, one may obtain the same result by assuming 
that the individual fiber centers follow the macroscopic 
strain field in an affine way. Then, relative displacements 
of centers of neighbouring fibers would be proportional 
to their typical distance. This is of the order filament 
length If and again Sz cclf. Note, however, that the as- 
sumption of affine displacement of the fiber centers can- 
not be literally true for fibers intersecting at very small 
angles 9^0. To avoid a diverging transverse deflection 
yi = —6z cot 6 —^ oo the two fibers will most likely not 
experience any relative motion at all and Sz 0. Truly 
affine displacements can therefore only be established on 
scales larger than the filament length. It should also be 
clear, that the assumption of afhne displacements of the 
fiber centers is different from the usual approach of as- 
signing afhne deformations on the scale of the single seg- 
ment. The latter would lead to deformations S^^e oc Is, 
proportional to the length Ig of the segment. Instead, 
axial displacements of the fiber as a whole are, by con- 
struction of the floppy mode, directly translated into non- 
afftne deformations Sna. cclf, which do not depend on the 
length of the segment. 

Restoring the radius r to its finite value, the floppy 
modes acquire energy and lead to bending of the fibers. 
A segment of length Is will then typically store the en- 
ergy Wb{ls) — ~ 1^1} /l^- By averaging over the 
segment length distribution Eq.lO one may calculate the 
average bending energy (Wb), stored in a fiber consisting 
of n ~ plf segments. 



{Wb) ^ plf 



dlsP{ls)^ 



(3) 



We assume the integral to be regularized by a lower cut- 
off length Zmin, that we now determine in a self-consistent 
manner. Physically, Imi-a corresponds to the shortest seg- 
ments along the fiber that contribute to the elastic en- 
ergy. Even though we know (see Eq.JSjl) that arbitrar- 
ily short segments do exist, their high bending stiffness 
k±(ls) oc Ij^ makes their deformation increasingly ex- 
pensive. Segments with length Is < Zmin will therefore 
be able to relax from their floppy mode deformation Sna, 
thereby reducing their bending energy from Wb^lmin) to 
nearly zero. However, due to the length constraints this 
relaxation necessarily leads to the movement of an en- 
tire secondary fiber and to the excitation of a floppy 
mode there. By balancing ■u;f,(/i„in) = (Wb) this gives 
Zniin — ^/p^lf and for the average bending energy of a 
single fiber (Wb) — K/lf{plf)^. This implies for the mod- 
ulus Gna — p/lf{Wb) oc p^ , which compares well with the 
simulation result of fi — 6.67. What is more, by equat- 
ing the energy (Wb) with (Wg) — Klfr~^ valid in the 
affine stretching regime, one can also infer the crossover 
exponent ly = 3. 

In summary, we have succeeded in explaining the elas- 
ticity of the bending dominated regime starting from the 



microscopic picture of the floppy modes that character- 
ize directly the deformation field deep inside the non- 
affine regime. Alternatively, one might try to understand 
the emergent non-affinity in a perturbative approach that 
considers deviations from an affine reference state. Such 
a line of reasoning has recently been suggested in [l8| . 
where non-affine boundary layers, growing from the fila- 
ment ends, are assumed to perturb the perfect affine or- 
der. However, comparing with their simulation data the 
authors could not confirm the scaling picture unambigu- 
ously and acknowledged the need for further numerical 
as well as improved theoretical work [l^ . Thus, non- 
affine elasticity in fibrous networks appears to be intrin- 
sically a non-perturbative strong-coupling phenomenon 
for which the floppy mode picture provides the correct 
low-energy excitations. As we will explicitly show next, 
one particular strength of our approach is that the scal- 
ing picture can readily be extended to a full theory that 
self-consistently calculates the modulus in a non-affine 
effective medium theory. 

To set up the theory we consider a single filament to- 
gether with its cross-links that provide the coupling to 
the medium. The energy of this assembly consists of two 
parts. First, the bending energy of the primary fiber 



«'.i»wi = f/(gl)'«^. 
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due to a transverse deflection y{z). A second "stretch- 
ing" energy contribution arises whenever a cross-link 
deflection yi = y{zi) differs from its prescribed value 
iji — —Sz cot 9i and may be written in the form of an 
harmonic confining potential Ws{yi) = \ki{yi — yiY that 
acts individually on each of the n plf cross-links. It 
allows the filament to reduce its own energy at the cost 
of deforming the elastic matrix it is imbedded into. Per- 
forming a configurational average (.) over cross-link posi- 
tions Zi and orientations 9i we obtain the average elastic 
energy stored in a single fiber as 



{W) ^ {myn(Wb[y{z)] 
\y(^-) \ 



" k 




(5) 



To solve the model we further need to specify the stiff- 
ness hi ~ k{9i) of the medium that relates to the relax- 
ation mode of a cross-link on the primary filament from 
its fioppy mode deflection. As we have argued above, any 
relaxation of this kind must act as axial displacement on 
a secondary fiber, thus exciting a new floppy mode there. 
The energy scale associated with this is (W) such that 
we can write 



k{9,) = 2{wy- 



5z^ 



(6) 



where the angular dependence derives from the projec- 
tion onto the axis of the secondary filament. Eqs.© and 
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FIG. 2: Graphical solution of Eqs.0 and Q for various num- 
bers n of cross-links obtained by calculating the intersection 
between the left side of the equation {W)ihs (bisecting line, 
dashed curve) with the right side (VF)rhs (full curves). The 
different curves for a given n correspond to ensembles of vary- 
ing size. They seem to diverge in the limit (W^)rhs5 ^ (W^)ihs. 
In fact, there (and only there) the averaging procedure is ill 
defined |2^ . Inset: Resulting dependence of (W) on n. 



© represent a closed set of equations to calculate the 
configurationally averaged deformation energy (W) as a 
function of the number of cross-links n. In implement- 
ing this scheme we have generated ensembles of filaments 
with a distribution of cross-linking angles as given by 
Eq.(^ and segment-lengths according to Eg.©. Note, 
that there is no free parameter in this calculation. The 
equations are solved graphically in Fig[21by plotting both 
sides of Eq. ((S]) as a function of {W) . The point of inter- 
section, which solves the equation, is shown in the inset 
as a function of the number of cross-links n. For the 
same parameter window as used in the network simula- 
tions ^6|, it yields the scaling behaviour of (W) oc n^ ''^. 
This implies for the modulus the exponent of ^ = 6.75, 
which improves upon the simple scaling picture presented 
above and provides a very accurate calculation of the 
scaling exponent fi. 

In conclusion, we have succeeded in deriving the 
macroscopic elasticity of random fibrous networks start- 
ing from a microscopic description of the displacement 
field in a manner that does not rely on the notion of affine 
deformations. We have given a floppy mode construction 
that may be applied to any two or three-dimensional net- 
work with fibrous architecture, for example paper or bi- 
ological networks of semiflexible filaments. It may also 
be shown to be relevant to systems where the constraint 
of straight fibers is relaxed The unusually strong 

density dependence of the modulus found here is a con- 
sequence of the exponential segment length distribution 
(PI and the presence of the length-scale Zmin ■ While iden- 
tification of the floppy modes has been recognized to be 
highly important for a description of force transmission 



in granular media or the jamming transition in colloidal 
systems, one can rarely give the exact form of these zero- 
energy excitations. On the contrary, we have achieved an 
explicit construction of the floppy modes that can be put 
in the form of localized elementary excitations ("floppi- 
ons" ) affecting only single filaments and their immediate 
surroundings. 

It is a pleasure to acknowledge fruitful discussions with 
David Nelson and Mikko Alava. 
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Due to the finite fiber length there are also two- and 

three-fold coordinated cross-links. 

A similar construction holds in 3d, where, in addition, 
each cross-link acquires a floppy "out-of-plane" degree of 
freedom. 



